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ff^ ■ Abstract 

O ■ We compute the upper bound on the mass of the hghtest Higgs boson in the 

Qs . Minimal Supersymmetric Standard Model in a model-independent way, includ- 

fH ■ ing leading (one-loop) and next-to-leading order (two-loop) radiative corrections. 

C^, We find that (contrary to some recent claims) the two-loop corrections are nega- 

Qh! tive with respect to the one-loop result and relatively small (< 3%). After defining 

(— I I physical (pole) top quark mass Mt, by including QCD self-energies, and physical 

Higgs mass Mh, by including the electroweak self-energies H (M^) — n(0), we 

obtain the upper limit on Mh as a function of supersymmetric parameters. We 

H ' include as supersymmetric parameters the scale of supersymmetry breaking Ms, 

the value of tan /? and the mixing between stops Xt = At + fi cot /3 (which is re- 
sponsible for the threshold correction on the Higgs quartic coupling) . Our results 
do not depend on further details of the supersymmetric model. In particular, for 
Ms < 1 TeV, maximal threshold effect Xf = 6MJ and any value of tan /?, we 
find Mh < 140 GeV for Mt < 190 GeV. In the particular scenario where the top 
is in its infrared fixed point we find Mh < 86 GeV for Mt = 170 GeV. 



CERN-TH.7334/94 
July 1994 



*Work partly supported by CICYT under contract AEN94-0928, and by the European Union under 
contract No. CHRX-CT92-0004. 

T Supported by a grant of Coniunidad de Madrid, Spain. 

*0n leave of absence from International School of Advanced Studies, ISAS, Trieste. 



1 Introduction 

There are good reasons to believe that the Standard Model (SM) is not the ultimate 
theory since it is unable to answer many fundamental questions. One of them, why 
and how the electroweak and the Planck scales are so hierarchically separated, has mo- 
tivated the proposal of the Minimal SupersjTumetric extension of the Standard Model 
(MSSM) as the underlying theory at scales of order 1 TeV. Indeed, the supersymmet- 
ric scale cannot be very large not to spoil the solution to the hierarchy problem, but 
the ambiguity about its maximal allowable size makes it very difficult to give model- 
independent predictions testable in accelerators. 

Fortunately, the supersymmetric predictions for the upper bound on the lightest 
Higgs boson mass represent a (perhaps unique) exception to this rule. Therefore, they 
are crucial for the experimental verification of supersymmetry. This importance is 
reinforced by the fact that uncovering the Higgs boson of the SM is one of the main 
challenges for present (LEP, Tevatron) and future (LEP-200, LHC) accelerators. Much 
work has been recently devoted to this subject [1-10], but still there is a substantial 
disagreement on the final results (see e.g. [5,9,10]), especially concerning the two-loop 
corrections (which, as we will see, are crucial). 

The aim of this paper is to evaluate the two- loop predictions of the MSSM on the 
lightest Higgs mass in a consistent and model- independent way. In fact, our results 
will not depend on the details of the MSSM {e.g. the assumption or not of universality 
of the soft breaking terms, gauge and Yukawa unification, etc.). Actually, when we will 
refer to the MSSM, we will simply mean the supersymmetric version of the SM with 
minimal particle content. The consistency of the approach will allow us to show up the 
reasons of the disagreement between previous results. We also give the predictions on 
the Higgs mass for a particularly appealing scenario, namely the assumption that the 
recently detected top quark [^ has the Yukawa coupling in its infrared fixed point. It 
will turn out that in this case the Higgs boson should be just around the corner. 

Let us briefiy review the status of the subject and the most recent contributions on 
it. The MSSM has an extended Higgs sector with two Higgs doublets with opposite 
hypercharges: Hi, responsible for the mass of the charged leptons and the down-type 
quarks, and H2, which gives a mass to the up-type quarks. After the Higgs mechanism 
there remain three physical scalars, two CP-even and one CP-odd Higgs bosons. In 
particular, the lightest CP-even Higgs boson mass satisfies the tree-level bound 

m^ < M| cos^ 2/3, (1.1) 

where tan/5 = ^2/^1 is the ratio of the Vacuum Expectation Values (VEV's) of the 
neutral components of the two Higgs fields if 2 and Hi, respectively. The relation 
(|1.1| ) implies that m"^ < M|, for any value of tan/3 which, in turn, implies that 
it should be found at LEP-200 [Q. However, the tree level relation (|1.1|) is spoiled 
by one-loop radiative corrections, which were computed by several groups using: the 
effective potential approach 0], diagrammatic methods [p[ and renormalization group 



(RG) techniques 0]. All methods found excellent agreement with each other and large 
radiative corrections, mainly controlled by the top Yukawa coupling, which could make 
the lightest CP-even Higgs boson to escape experimental detection at LEP-200. In 
particular, the RG approach (which will be followed in the present paper) is based on 
the fact that supersymmetry decouples and below the scale of supersymmetry breaking 
Ms the effective theory is the SM. Assuming M| -C M'^ the tree-level bound ( |1.1D is 
saturated at the scale Ms and the effective SM at scales between Mz and Ms contains 
the Higgs doublet 

H = Hicosp + ia2H;smP, (1.2) 

with a quartic coupling taking, at the scale Ms, the (tree level) value 

X = \ig' + g")cos'2p. (1.3) 

In these analyses [0] the Higgs mass was considered at the tree-level, improved by one- 
loop renormalization group equations (RGE) in the 7- and /5-functions, thus collecting 
all leading logarithm corrections. 

Since the relative size of one-loop corrections to the Higgs mass is large (mainly for 
large top quark mass and/or small tree level Higgs mass) it was compelling to analyze 
them at the two-loop level. A first step in that direction was given in ref. [^ where 
two-loop RGE improved tree level Higgs masses were considered. It was found that 
two-loop corrections were negative and small. In ref. the Higgs mass received all 
leading logarithm and part of the next-to-leading logarithm corrections. Subsequent 
studies of the effective potential improved by the RGE p, |^, H have shown that the 
L-loop improved effective potential with (L-l-l)-loop RGE is exact up to Lth-to-leading 
logarithm order |j^Q. This means that for fully taking into account all next-to-leading 
logarithm corrections the one-loop effective potential (improved by two-loop RGE) is 
needed. 



Finally, two papers 0, [T^ have recently appeared aiming to refine the two-loop 
analysis of ref. . The authors of ref. ^ , following the RG approach, find that the 
two-loop correction to the Higgs mass is positive and sizeable, whereas the authors 



of ref. |T^, following diagrammatic [^ and effective potential approaches in the 
framework of the MSSM with various approximations, find that the two-loop correction 
to the Higgs mass is negative and sizeable in contradiction with ref. [§. As we will see 
in detail in section 5 we are in disagreement with the results of ref. ||^ (showing how a 
correct treatment of all the relevant effects would make the results presented in ref. 
to agree with ours), and in agreement with the overall result of ref. |]10[ (though not 
with the relative size of the two-loop corrections). 

In this paper we use the RG approach and the SM one-loop improved effective 
potential with two-loop RGE to compute the lightest Higgs mass of the MSSM up to 
next-to-leading logarithm order. At this level of approximation, one has to keep control 



"'^Strictly speaking this has been proven M for a theory with a single mass scale. 



on the validity of the perturbative expansion for V. In particular, although the whole 
effective potential is scale- invariant, the one- loop approximation is not, thus one has to 
be careful about the choice of the renormalization scale. As we will see this fact is at 
the origin of some misunderstandings in previous works. Furthermore, since the Higgs 
mass is computed to one-loop order there are a number of one-loop effects that need 
to be considered for the consistency of the procedure: 

i) As it is explained in Appendix A, the tree level quartic coupling (pT3|) receives 
one-loop threshold contributions at the Ms scale. These are given by 

where ht is the top Yukawa coupling in the SM and 

Xt = At + ficotp (1.5) 

is the stop mixing. 

The correction (|1.4| ) has a maximum for X^ = 6M|. For that reason, in our 
numerical applications we will take two cases: Xt = 0, i.e. no mixing, and X^ = 6M|, 
i.e. maximal threshold effect. Notice also that X^ = 6M| is barely consistent with the 
bound from color conserving minimum [^], so the case of maximal threshold really 
represent a particularly extreme situation. 

In addition to the previous effect, there appear effective higher order operators 
(-D > 6), which for Ms > 1 TeV turn out to be negligible (for details see Appendix A). 

ii) One-loop contributions to the top quark self-energy relating the running top 
mass rrit to the (physical) propagator pole top mass Mj. We will find this effect gives 
a negative and sizeable contribution to the Higgs mass 

iii) One-loop contributions to the Higgs self-energy relating also the running Higgs 
mass rriH to the physical Higgs mass M^ (for details see Appendix B). We will find 
scale dependent contributions to the self-energy, removing the scale dependence of Mfj, 
and other scale independent contributions. This effect will be found to be positive in 
all cases, partially cancelling the previous effect. 

The contents of this paper are as follows: In section 2 we present and analyze the 
one-loop effective potential, for the effective SM in the range between Mz and Ms, 
improved by the RGE. We compare different treatments of the effective potential and 
define the running Higgs masses. Our treatment of the effective potential is presented 
in detail in section 3, and the numerical results are presented in section 4, where we 
make use of effects (i)-(iii) above. Numerical and conceptual comparison with other 
recent approaches is presented in section 5, and our conclusions are drawn in section 
6. In Appendix B we present analytic expressions for the one- loop contribution to the 
Higgs boson self-energy H(g^) — H(0), computed in the 't Hooft-Landau gauge and 
using the MS renormalization scheme. We also comment on the introduction of the 
Higgs sector both in the effective potential and in the Higgs boson self-energy. 



2 The one-loop effective potential 



Our starting point in this section will be the effective potential of the Standard Model. 
This can be written in the 't Hooft-Landau gauge as f^ 



v,ff{fi{t),\,{ty,4>{t)) = Vo + Vi + 



(2.1) 



where Aj = {g,g',X,ht, m^) runs over all dimensionless and dimensionful couplings and 
Vo, Vi are respectively the tree level potential and the one-loop correction, namely 
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where we have used the MS renormalization scheme |jT2|. The parameters A(t) and 
m{t) are the Standard Model quartic coupling and mass, running with the RGE, while 
the running Higgs field is 

m = m<pc. (2.4) 



being the classical field and 



e(t)=e-/o^(*')'^^ 



(2.5) 



where 7(t) is the Higgs field anomalous dimension. Vi in eq.( |2.3D contains the radiative 
corrections where only the top quark and the W, Z gauge bosons are propagating in 
the loop. For the moment we will disregard those coming from the Higgs and Goldstone 
boson propagation]^. They can be easily introduced, as shown in Appendix B, without 
altering the numerical results that are obtained in this paper. The mass parameters in 
(|27^) are given by 



m|(t) 
ml{t) 
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(2.6) 
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him^t), 



^In particular, the Goldstone boson contribution to the effective potential generates, on the running 
Higgs mass, an infrared logarithmic divergence. This divergence is cancelled when the physical Higgs 
mass is considered, as will be shown in Appendix B. 



where g, g' and ht are the S'f/(2), f/(l) and top Yukawa couphng, respectively. Finally 
the scale /i(t) is related to the running parameter t by 

^^{t) = fie\ (2.7) 

where /x is a fixed scale, that we will take equal to the physical Z mass, Mz- 

The complete effective potential Veff{fi{t), Aj(t); (pit)) and its n — th derivative are 
scale-independent (see e.g. ref. P), i.e. 

^^0, (2.8) 

where 

Pin 

Ve?} ^ C{t)^—-V^fMt), A.(t); 0(t)). (2.9) 



The above property allows in principle to fix a different scale for each value of the 
classical field (pc, i-e- fJ'it) = /(0c) or, equivalently, t = t{(j)c) = log(/(<^c)//^)- In that 
case Veff becomes Veff{4>c) = Veff{f{(pc), Aj(t(0c))- Then, from ( WM and ( p. 91) one can 
readily prove that 



eff 



= '^^^, (2.10) 



as expected. This procedure has been used in some previous worksQ. However, one 
should notice that even though the whole effective potential is scale invariant, the one- 
loop approximation is not. Therefore one would need a criterium to fix the function 
f{(pc)' the only possible one would be to minimize the radiative corrections improving, 
so, perturbative expansion. Assuming that only one field, say / with squared mass 
■m'j{t), is contributing to the one- loop radiative corrections, or that the latter are dom- 
inated by this field, the most natural choice would be /i^(t) = m'j{t). Nevertheless this 
approach has (in our context) some drawbacks: 

i) The t-dependence oim'j{t) is implicit through the RGE and so the function t = t{(j)c) 
cannot be explicitly written. 



ii) In most cases, in particular in our one- loop correction ( |2.3| ), one cannot assume, 
within the required degree of accuracy, that loop corrections are dominated by 
only one field. 

r(n). 



iii) When computing V^^l{(j)c) in (|2.10| ) the vacuum energy Q{Xi{t) , fi{t)) has to be 



specified since it also acquires a 0c-dependence through the variable change 
t = t{(j)c)- Furthermore, the whole dependence on 0c becomes now much more 
involved. 



^E.g. in studies of the effective potential stability B to control large logarithms that can arise for 
large field values. It has also been used in ref. j^, but as we will see, it is an unnatural complication 
in this context. 



iv) The last, but not the least, drawback is that once one has chosen a particular 
function t(0c), one looses track of scale invariance and so there is no way of 
checking how good the approximation is at the minimum of the effective potential. 

For the above mentioned reasons we will keep t and 0c as independent variables. We 



will minimize the effective potential (2J.), truncated at one-loop, at some fixed scale 
t = t*, i.e. 

dVeff 



d(j){t* 



= 0, (2.11) 

0(t*)=Wt*)> 

which will determine the VEV (0(t*)). Our criterium to fix the scale t* will be explained 
in the next section. The scale independence of the whole effective potential implies that 



d (0(t)) 
dt ^(t) 



0. 



(2.12) 



Therefore, assuming that t* lies in the region where the one-loop approximation to 
the effective potential is reliable, the VEV of the field at any scale can be obtained 
throughfl 



m)) = m*)) 



at*y 



Accordingly, we must impose 



t;=(0(t^)) = (0(r))|M, 



where tz is defined as /i(t^) = Mz and 



246.22 GeV, 



(2.13) 



(2.14) 



(2.15) 



is the "measured" VEV for the Higgs field[| |T^ . 

We can trade {(pit*)) by m?(t*) from the condition of minimum ( 2.11|) , which trans- 
lates, using ( p.l4| ), into the boundary condition for rn?{t): 
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(2.16) 



2/i2(t*) 

^Notice that (<j){t)) defined by ( 2.13|) coincides with the value of 0(t) which minimizes the whole 
effective potential. In the one-loop approximation that we are using this is no longer true, as we will 
see later on. 

^We are neglecting here one-loop electroweak radiative corrections to the muon /3-decay slightly 
modifying the relation ( 2.15 ), see e.g. ref. ||l5|. We thank S. Peris for a discussion on this point. 



We can define now innit*) as the second derivative of tlie effective potential at tlie 
minimum, evaluated at the scale t*, i.e. 
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(2.17) 



where we have used ( |2.16| ). At an arbitrary scale t we can use the scale independence 
property of the whole potential, see (|2^ ), and write 



(2.18) 



m|(t)=m^(r)^. 



which will be our definition of running mass. 

In the region where our approximated effective potential is scale invariant, definition 



(|2.18| ) should be equivalent to taking the second derivative with respect to (pit), evaluate 
it at {(l){t)) and use, for the VEV, the relation ( |2.13|) , i.e. 
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(2.19) 



Of course, exact scale invariance would imply that ( |2.18| ) and (|2.19|) are equal. This 
will allow to cross check the reliability of our approach. 



3 Our approach 

In this section we will describe our approach to the problem of determination of the 
Higgs mass (|2.18 ). The effective potential is written in (|2.2|) and ( p. 3D , and the pa- 
rameters on which it depends, A(t), g(t), g'(t), g^it), ht(t), ^(t) and im?{t), satisfy a 
system of coupled RGE with t-dependence governed by /3a, (3g, Pg', Pg^, Pht: 7; Pm'^y 
i.e. the corresponding /5-functions and anomalous dimension of the Higgs field, which 
are evaluated to one- or two-loop order |§]. In fact, we will often consider two cases: 



a) The one-loop case, where (3 and 7-functions are considered to one-loop and the 

effective potential is approximated by the tree level term (|2.2|) . In this case the 
leading logarithms are resummed to all-loop in the effective potential. The Higgs 
mass includes then all-loop leading logarithm contributions |^, 0]. 

b) The two-loop case, where jS and 7-functions in the RGE are considered to two-loop 

order and the effective potential is considered in the one-loop approximation, ( p.2|) 
and (|2.3| ). In this case the leading and next-to-leading logarithms are resummed 
to all-loop in the effective potential. The Higgs mass includes then all-loop leading 
and next-to-leading logarithm contributions [0. 

We will impose as boundary conditions: 

• For the gauge couplings: 

g{Mz) = 0.650, 

g'{Mz) = 0.355, (3.1) 

gs{Mz) = 1.23 . 



For the top Yukawa coupling: 



h,{m,) = ^^ , (3.2) 

V 



where rrit is the running mass: ■mt{^(t) = rrit) = rrit. 
For the quartic coupling: 



1 



4 



. ^ , 3ht{Ms) fX^ X, 



X{Ms) = - g\Ms) + g'\Ms) cos^ 2(3 + ;; f' 2-^ - — ^ . (3.3) 



-4 



167r2 \ Mi 6M; 



• For the mass m?{t) we will take as boundary condition the value determined by 
eq. (|2.16|) at the scale t* . 

Note that the previous boundary conditions depend on the values of m^, t* {i.e. the 
minimization scale) and the supersymmetric parameters Ms, Xt and tan/3. 

Our main task now will be to determine the optimal minimization scale t* as a 
function of mt. Our criterion will be that t* be in the region where the effective 
potential is more scale independent. We estimate this in the following way. For fixed 
t* and mj, all boundary conditions ( |3.1|) -( p73D are determined and the effective potential 
(121] ) can be computed for any value of t. Thus the corresponding value of 0(t), say 
0min(^); that minimizes the one-loop potential can be numerically evaluated. Had we 
considered the whole effective potential, its scale independence property (|2.12| ) would 



imply that (pminit) / ^{t) = {(f){t))/C,{t), where (0(t)) has been defined in ( |2.13| ), is t- 
independent. Since the one- loop approximation is not exactly scale independent the 
most appropriate scale is the scale tg for which 0min(^)/^(^) has a stationary point, i.e. 



d (j)mm{t) 



= 0. (3.4) 

t=ts 



dt i(t) 

Of course the scale at which (|3.4D occurs depends on t* and m^, i.e. tg = ts{t*,mt 



Therefore, for a given value of rrit the optimal scale t* is defined as 

ts = ts{t*,mt)=f, (3.5) 



i.e. it is the scale which simultaneously defines the boundary condition ( p.l6|) and 
agrees with the extremal of the function 0min(^)/^(^)- 

This procedure is illustrated in Fig. 1, where we plot 4>min{t)/C,(t) vs. /i(t) for 
rrit = 160 GeV and supersymmetric parameters Ms = 1 TeV, Xt = 0, cos^2j3 = 1. 
The solid curve shows the two-loop result whose stationary point satisfies eq. ( p.5|) 
and therefore defines the optimal scale t*. Note that there is no fine tuning in the 
choice ( |3.5| ) since any point near the stationary point would be equally appropriate 
because the curve (f)mm{t) / ^{t) is very fiat in that region. We have also shown in Fig. 1 
the curve (f)min{t)/^{t) in the one-loop approximation (dashed curve) with fi{t*) fixed 
by the two- loop result. We can see that the one- loop curve is much steeper than the 
two-loop curve, which means that the one-loop result, i.e. the tree level potential ( [^.2[ ) 
improved by the one-loop RGE, is far from being scale independent at any scale. This 



feature has also been observed for the MSSM one-loop effective potential [|T7|. 



We plot in Figs. 2a,b,c,d, /u(t*) as a function of rrit for the different values of su- 
persymmetric parameters Ms = 1, 10 TeV; X^ = 0, 6M|. In all plots the solid curve 
corresponds to cos^ 2/5 = 1 and the dashed curve to cos^ 2/3 = 0. 

Once we have determined n{t*) for fixed values of rrit and all supersymmetric pa- 
rameters, the Higgs running mass m//(t) is given by ( p^.lS] ), while the mass defined as 



the second derivative of the effective potential rnH,der{t) is given by (|2.19| ). By definition 
both masses coincide at t* 

mnit*) =mH,der(t*), (3.6) 

and the ratio mH4er(t)/mH(t) should be equal to one for an exactly scale independent 
effective potential. Consistency of our procedure requires the curve 'rnH,der{t)/mHit) to 
be flat in the region where we minimize, i.e. at t*. In Fig. 3 we plot rriH^derif) / rriHit) 
for m^ = 160 GeV and the values of supersymmetric parameters as in Fig. 1, Ms = 1 
TeV, Xt = 0, cos^ 2/5 = 1. We see that at the point /i(t*) the curve is in its flat region, 
though /i(t*) does not exactly coincide with the extremum oi mH,der(t)/mH(t). 

Neither mn nor rrit are physical masses. They are computed from the effective po- 
tential, i.e. at zero external momentum, and need to be corrected by the corresponding 
polarizations to obtain the propagator pole physical masses. This will be done in sec- 
tion 4. For the time being, and just to compare with other results in the literature, we 



will neglect the shift to the physical poles and define the top mass by the boundary 



condition ( |3.2|) and the Higgs mass by the usual condition 

niH {fi(t) = uih) = niH. (3.7) 

In this way the Higgs mass tuh can be unambiguously determined. We plot in Fig. 4 
rriH as a function of rrit for the values of supersymmetric parameters Ms = 1 TeV, 
Xt = 0, cos^ 2/5 = 1. The thin solid line corresponds to the two-loop result and the 
dashed-line corresponds to the one- loop result. (We will disregard for the moment the 
thick solid line, which corresponds to shifting the running masses to propagator poles.) 
The main feature that arises from Fig. 4 is that the two-loop corrections are negative 
with respect to the one- loop result. This feature is in qualitative agreement with our 
previous two-loop result (where the one-loop corrections to the effective potential 
where not considered) and with others from different authors [0. This comparison 
will be done in some detail in section 5. 

4 Numerical Results 

We will present, in this section, the numerical results on the Higgs mass, evaluated 
in the next-to-leading logarithm approximation, as a function of the top quark mass. 
The running top quark mass rrit that we have been using in the previous section was 
evaluated at the scale /i(t) = rrit, i-^- fntifnt) = iT^t (see eq. (|3.2|) ). However the running 
mass does not coincide with the gauge invariant pole of the top quark propagator M^. 
In the Landau gauge the relationship between the running rut and the physical (pole) 
mass Mt is given by 



M, 



■^^4«,(M,) 



3 n 



mt{Mt). (4.1) 



On the other hand, the running Higgs mass, mnit), given by eq. ( |2.18|) , has a scale 
variation 

^^ = 27<(t). (4.2) 

The propagator pole Mh is related to the running mass through (see Appendix B) 

Ml = m%{t) + ReH(M|) - ReH(O), (4.3) 

where H(g^) is the (renormalized) self-energy of the Higgs boson. In our calculation 
(case (b) in section 3) it is enough to compute the Higgs self-energies at the one-loop 
level since we are computing the Higgs masses to one-loop. The imaginary part of 
H(M|f) — H(0) contributes to the Higgs width. Assuming Mh < 2Mw, the Higgs is 
stable at tree level, and eq. (|4.3|) reads 



M^ = m| + H(M^)-H(0). (4.4) 

10 



The calculation of n(M^) — 11(0) is presented in Appendix B. We note that the scale 
dependence (at the one- loop level) of m//(t), as provided by ( |4.2| ), is cancelled by the 
scale dependence of n(M^) — 11(0). The only remaining scale dependence of M^ comes 
from two-loop contributions. This feature is exhibited in Fig. 5 where we plot Mh in 
the range of /i(t) between Ms and Mz for rrit = 160 GeV and the supersymmetric 
parameters Ms = 1 TeV, Xt = 0, cos^2/3 = 1. We can see that Mh has a negligible 
variation (~ 0.5 GeV) between Ms and Mz- For the sake of comparison we have also 
plotted the running mass mH{t) whose variation is more appreciable (~ 2.5 GeV). This 
effect is more accentuated for larger top masses. 

We can see from ( [4.1D and from Fig. 4 that the effect of considering the pole mass Mt 
is negative on the two-loop corrections to the Higgs mass, while the effect of n(M^) — 
n(0) is positive (~ 2 GeV for mt = 160 GeV and ~ 5 GeV for rrit = 215 GeV) thus 
partially compensating the effect of ([4.1|) . The global effect is negative and small as 
can be seen in Fig. 4. The thick solid line indicates Mh as a function of Mf. It is below 
the thin solid line which was the two- loop evaluation of ttih as a function of rrit. The 
comparison with the one-loop result (dashed line) shows that two-loop corrections are 
negative with respect to the one- loop result. Numerically they are small (~ 1 GeV) 
for Mt = 120 GeV and larger (~ 6 - 7 GeV) for Mt = 220 GeV. 

In Fig. 6a,b,c,d we plot Mh as a function of Mt for values of supersymmetric 
parameters Ms = 1, 10 TeV; X"^ = 0, 6M|; cos^ 2/3 = 0, 1. In all cases the solid curve 
corresponds to cos^ 2/? = 1 and the dashed curve to cos^ 2/3 = 0. Notice that the 
dependence on the mixing parameter Xt is sizeable. In all the figures we have used the 
lower limit on the top quark mass, Mt > 120 GeV at 95% C.L., from the CDF dilepton 



channel [^ . If we use the recent evidence for the top quark production at CDF with 
a mass Mt = 174 ± 101^2 GeV ^^ and the bounds for supersymmetric parameters 



Ms < 1 TeV and maximal threshold correction X"^ = 6M|, we obtain the absolute 
upper bound Mh < 140 GeV. 

The dependence of Mh on tan/3 is exhibited in Fig. 7 where we fix Mt = 170 GeV 
and Ms = 1 TeV. The solid curve corresponds to the absolute upper bound for the 
mixing Xt = 6M| and the dashed curve to the case of zero mixing. Concerning the 
dependence of Mh on the stop mixing Xt parameter (or equivalently the stop splitting), 
we have found it to be sizeable, as can be seen from Figs. 6. 

5 Connection with other approaches 

Some papers have recently appeared aiming to estimate the mass of the lightest Higgs 
boson up to next-to-leading order in the MSSM, and with apparently contradictory 
results. In this section we will comment on those papers in the context of our formalism 
and will exhibit the origin of the disagreements. 

In ref. the RG approach was used to estimate the mass of the lightest Higgs boson up 
to next-to-leading order. The relevant points of their calculation are the following: i) 
They considered the one-loop correction to the effective potential in the approximation 
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g = g' = 0, ii) They neglected the wave- function renormahzation leading to physical 
masses for the top quark and Higgs boson, iii) They minimized the effective potential 
at the scale fi{t*) = v. In addition they only consider the case of zero stop mixing, 

Contrary to our results, they found that the two-loop correction is positive and 
sizeable with respect to the one- loop result. For instance in the typical case Ms = 1 
TeV, Xt = 0,cos^2/? = 1, they find for rrit = 170 GeV the two-loop result rriH ~ 117 
GeV, while we find for Mj = 170 GeV the two- loop result Mh ~ HI GeV. We have 
been able to trace back the difference between the two results to the points i)-iii) 
above. In particular, the authors of ref. [Q find for the same values of the parameters 
a positive two-loop correction with respect to the one-loop result ~ 7 GeV while our 
two-loop result is smaller than the one-loop result by ~ 3 GeV. To understand the 
origin of this discrepancy we have plotted in Fig. 8 the physical Higgs mass Mh as 
a function of the minimization scale /i(t*) for Mt = 170 GeV, Ms = 1 TeV and 
X^ = 0. Thick lines correspond to cos^ 2/3 = 1 and thin lines to cos^ 2/5 = 0. Solid 
lines represent M^ evaluated in the two-loop approximation and dashed lines in the 
one- loop approximation. Our chosen value of /i(t*) is indicated in the figure with an 
open diamond and that chosen in ref. 0, fi(t*) = v, with an open square. We can 
see that Mh evaluated at two-loop is very stable against /i(t*): it varies ~ 5 GeV 
in the whole interval Mz < /^(t*) < Ms and ~ 2 GeV between the diamond and 
the square. However the value of Mh evaluated at one-loop is very unstable. In fact 
the two-loop correction changes from negative at /i(t) = Mz (^ 10%) to positive at 
fi{t) = Ms (> 20%). In the region of our chosen value of A*(t*) it is negative and ~ 3 
GeV while at yu(t*) = f it is positive and greater ~ 7 GeV. Fig. 8 shows that though 
our choice of fi{t*) is more in agreement with perturbation theory than fi{t*) = v the 
difference is however not important when plotting the physical Higgs mass (which was 
not considered in ref.0). 



In ref. |T0| the diagrammatic and effective potential approaches were used to evaluate 
the lightest Higgs mass at two-loop order in the framework of the MSSM. Various 
approximations, like g = g' = in the effective potential, were used, and only the 
case of zero stop mixing, X^ = 0, was considered. Our results agree with those of 
ref. fig within less than ~ 3 GeV. In fact for Mt = 150 GeV, Xt = 0, cos^ 2/? = 1 



and Ms = 1 TeV {Ms = 10 TeV) ref. H finds Mh ~ 107 GeV {Mh ~ HO GeV) 



while we find from Figs. 6a and 6b Mh ~ 103 GeV {Mh ~ 110 GeV). We consider 
this agreement as satisfactory given the approximations used in ref. |jlO|. The fact that 



two-loop corrections found in ref. |T^ are sizeable with respect to the one-loop result 
can be explained as a consequence of the approximation used there to estimate the one- 
loop result. Taking Ms = 10 TeV and Mt = 150 GeV, they found mn = 138 GeV at 
one-loop using a simple approximation that takes into account only the leading part (~ 
M^ log{Mg/M^)) of the corrections, as it is common practice in some phenomenological 
analysis. Actually, a slightly more sofisticated approximation (as the one labelled 1(3 in 
their paper) or the one-loop result obtained by a numerical integration of the RGE (as 
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in our approach) gives mn ~ 115 GeV (both methods agree up to a difference ~ 2 — 3 
GeV due to subleading effects, such as those provided by considering gauge couphngs 
in the effective potential) for the above mentioned values of the parameters. This has 
to be compared with the two- loop result M^ ~ 110 GeV and shows that the net effect 
of two-loop corrections is indeed negative and smallR. 



Finally bounds on the lightest Higgs mass have been recently analyzed in ref. ||T^ 
in the context of models with Yukawa unification. The results of ref. |18[ are presented 
in the one-loop approximation. They choose Mz as the minimization scale and find 
larger bounds than in previous estimates. For instance it is found, for Mf = 170 
GeV and any value of the supersymmetric parameters such that Ms < 1 TeV, that 
Mh ^ 102 GeV. We have tried to reproduce their results using our formalism. In fact, 
their solution for small tan/3 is close to the fixed point solution where tan/3 and rrit 
are related through [|n 



mt = (196 GeV)[l + 2{a3{Mz) - 0.12)] sin/3. (5.1) 

Now, fixing Ms = 1 TeV and maximal mixing^ X"^ = 6M| we find, for Mt = 170 GeV 
that tan/3 = 1.74 from ( |5.1| ) and Mh ~ 99 GeV if we fix /i(t*) = Mz, in agreement 
with the result of ref. [|l8l. However this large value of Mh is a clear consequence 
of the chosen minimization scale where two-loop corrections are very large. In fact 
for /i(t*) = Mz our two- loop result gives Mh ~ 85 GeV. Moving to the region where 
perturbation theory is more reliable, as we have done along this paper, we would obtain 
the final two-loop result, Mh ~ 86 GeV. We have plotted in Fig. 9 Mh as a function 
of Mt for Ms = 1 TeV and tan/?, determined from ( ^.1|) , corresponding to the fixed 
point solution. The solid (dashed) line corresponds to the case of maximal mixing, 
Xt = 6M| (zero mixing, Xt = 0). 

6 Conclusions 

We have computed in this paper the upper bounds on the mass of the lightest Higgs 
boson in the MSSM including leading and next-to-leading logarithm radiative correc- 
tions. We have used an RG approach by means of a careful treatment of the effective 
potential in the SM, assuming that supersymmetry is decoupled from the SM, which 
we have shown to be an excellent approximation for Ms ^ 1 TeV or even much less. 
Our results have covered the whole parameter space of supersymmetric parameters; in 



^Notice that the results of ref. [Q only contain one- and two-loop leading and next-to-leading 
contributions to the Higgs mass while, as we have noticed before, our calculation of the Higgs mass 
includes leading and next-to-leading contributions to all-loop. This effect can be important for the 
case of large logarithms. We thank R. Hempfling for a discussion on this point. 

''An interesting result that can be confirmed from ref. pSl (see Fig. 10 in that paper) is that the 
maximal threshold limit is compatible with the non-existence of color breaking minima, as we already 
noted. For that reason removing the color breaking minima of the distribution plotted in Fig. 7 does 
not change the upper bound on the Higgs mass. 
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fact, when we refer to the MSSM, we simply mean the supersymmetric version of the 
SM with minimal particle content. We have also included QCD radiative corrections 
to the top quark mass, which gives a negative contribution to the Higgs mass as a 
function of the physical (pole) top mass Mf, and electroweak radiative corrections to 
the Higgs mass Mh, n(M^) — 11(0), providing a positive contribution to M^ and par- 
tially cancelling the former ones. Therefore the balance of including the contribution of 
QCD and electroweak self-energies to the quark and Higgs masses is negative and small 
(^ 2% for Mt < 200 GeV). We also made a reliable estimate of the total uncertainty 
in the final results, which turns out to be quite small (~ 2 GeV). 

Concerning the numerical results, we have found in particular that for Mt < 190 
GeV, Ms < 1 TeV and any value of tan/? we have M^ < 120 GeV for Xf = 0, and 
Mh < 140 GeV for X^ = 6M| (maximal threshold effect). We have also applied our 
calculation to a particularly appealing scenario, namely the assumption that the top 
has the Yukawa coupling in its infrared fixed point. Then, the bounds are much more 
stringent. E.g. for Mt = 170 GeV, Ms < I TeV and X^ = QM^, we find Mh < 86 
GeV. 

Two papers have recently tried to incorporate radiative corrections to the Higgs 
mass up to the next-to-leading order and with qualitatively different results. In ref. 0, 
using the RG approach, positive, and large next-to-leading corrections with respect to 



the one-loop results were found. In ref. |10], using diagrammatic and effective potential 



methods in a particular MSSM as well as various approximations, it was found that 
two-loop corrections are also sizeable, but negative with respect to the one-loop result! 

Using the RG approach, as ref. [P], we have found (unlike in ref. [^) that two- loop 
corrections are negative with respect to the one-loop result. We have traced back the 
origin of this disagreement in their choice of the minimization scale. Furthermore the 
authors of ref. ^ neglected various effects (as the contribution of gauge bosons to the 
one-loop effective potential, or the wave function renormalization of top quark and 
Higgs boson) and considered only the case with zero stop mixing. 

On the other hand, we have found that the abnormal size of the two-loop corrections 
obtained in [^ is a consequence of an excesively rough estimate of the one- loop result, 
but we are in agreement with their final two-loop result. In fact our two- loop results 
differ from those of ref. |]TU1 by less than 3%. Also our results show a large sensitivity 



of the Higgs mass to the stop mixing parameter. 

Finally we would like to comment briefly on the generality of our results. As 
was already stated, we are assuming average squark masses Mg ^ M^, and that all 
supersymmetric particle masses are > Ms. If we relax the last assumption, i.e. if some 
supersymmetric particles were much lighter, the value of the quartic coupling at Ms 



(see eq. ( p..3|) ) would be slightly increased and, correspondingly, our bounds would 
be slightly relaxed. We have made an estimate of this effect. Assuming an extreme 
case where all gauginos, higgsinos and sleptons have masses ~ Mz, we have found for 
Ms = 1 TeV and cos^ 2/5 = 1 an increase of the Higgs mass ~ 2%. For values of tan/? 
close to one (as those appearing in infrared fixed point scenarios) the corresponding 
effect is negligible. On the other hand, our numerical results have been computed 
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for Ms > 1 TeV. For values of Ms < 1 TeV the bounds on the hghtest Higgs mass 
are lowered. Hence, in this sense, all our results can be considered as absolute upper 
bounds. 
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Appendix A 



We comment here on the origin and size of threshold contributions to the quartic 
coupling of the SM, A, at the supersymmetric scale. We also evaluate the effective 
D=6 operators that are relevant for the Higgs potential, showing that for Ms > 1 TeV 
they are negligible. 

As it is well known, the one-loop correction to the MSSM effective potential is 
dominated by the stop contribution 

.2 



rMSSM 



^ 327r2 I 1 






rrin 



, mi 3 



(A.l) 



where mi, m,2 are the two eigenvalues of the stop mass matrix. In a good approxima- 
tion: 

m^ 2 = Ml + hlH2 ± htH2Xt (A.2) 

where Mg is the soft mass of the stops (we are assuming here that this is the same for 
the left and the right stops, which is a correct approximation since the RGE of both 
masses are dominated by the same QCD term); ht is the top Yukawa coupling in the 
MSSM and Xt = At + n cot (5 gives the mixing between stops {At is the coefficient of 
the soft trilinear scalar coupling involving stops and /i is the one of the usual bilinear 
Higgs term in the superpotential). Notice also that Ms is basically the average of the 
two stop masses, which is a reasonable choice for the supersymmetric scale when one 
is dealing with the Higgs potential, as it is our case. Thus we identify the scale Q at 
which we are evaluating the threshold effects with Ms- 

Now it is straightforward to obtain from ( |A.1|) the contribution to the H^ operator 
(recall that H is the SM Higgs doublet given by eq.(0|)). That is 



i{(f)^'-(ai)^'h'- <"' 

where have redefined ht to be the usual top Yukawa coupling in the SM. This gives a 
threshold contribution to the SM quartic coupling 
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Alternatively, eq. ( |A.4| ) can be obtained diagrammatically from two kinds of dia- 
grams exchanging stops (see the first paper of ref. 0|). Both methods give the same 
result. 

It is worth noticing that ( A.4 ) has a maximum for X'^ = 6M|, which means that 



the threshold effect cannot be arbitrarily large. Note also that X"^ = 6M| is barely 



consistent with the bound from color conserving minimum [^ , so the case of maximal 
threshold represents an extreme situation. 

Analogously, we can obtain from ( |A.1| ) the relevant effective D=6 operators, i.e. 



those proportional to H . These turn out to be 

which could also be obtained diagrammatically from four kinds of diagrams exchanging 
stops. It is easy to see that (|A.5|) produces negligible modifications in the process of 



electroweak breaking and in the Higgs mass. For example, for Ms = 1 TeV and 
maximal threshold in ( [A.4| ), i.e. X'f = 6M|, ( |A.5| ) gives modifications in the Higgs 



mass suppressed by a factor ~ 1/150 with respect to those induced by ( |A.4D . This also 



means that, regarding the Higgs potential, it is safe to decouple the MSSM from the 
SM for Ms = 1 TeV or even substantially smaller. 

Appendix B 

In this Appendix we shall discuss in more detail the relation between the running mass 
of the Higgs boson, extracted from the effective potential, and the physical Higgs mass 
and give the complete expression for the latter. 

The physical mass of the Higgs boson field M^ is defined as the pole of the prop- 
agator and it is both renormalization scheme and gauge (if calculated at all orders of 
perturbation) independent. 

We start with the Lagrangian 

= lid<l)Rf -^ml<l>l + ^6ZHid<l)n)' -l5m'<l>% + ... , (B.l) 

where 0o and 0/j = ZJj ' (f)Q = (1 + 5Zh)~^^'^4>o are the bare and renormalized Higgs 
fields, and ttiq and m^ are the bare and renormalized masses. With this convention 
?t1q = ?7i|, + 5m? — SZhttlJ^. Denoting by Fq (F/j) the inverse of the one-loop corrected 
bare (renormalized) propagator, we have 

Fi?,(p^) = ZhTo{p^) = Zh P^ -m%- 6m^ + 6ZHm\ - Ho(p^ 

= p' + 5ZHv'-{rnl + 5m^)-Yl^{p') , (B.2) 
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where no(p^) is the (unrenormahzed) self-energy of the Higgs boson field and in the 
last equality we have neglected higher order terms in the perturbative expansion. The 
renormalized self-energy is defined as 

liRip") = no(/) + 6m^ - dZup". (B.3) 

Thus 

Tn{p')=v'-[ml + Tln{p')) • (B.4) 

Consequently, the physical (pole) mass, M'jj, satisfies the relation 

Ml = ml + IiR{p^ = Ml) . (B.5) 

On the other hand, the running mass rnjj, defined as the second derivative of the 
renormalized effective potential (see eqs. ( |2.17| ) and ( [^.181 )), is given by 



d'K 



m 



eff 



H 



-rii(p' = 0)=m^ + n^,(p2 = 0) . (B.6) 



Comparing ( p.5|) and ( p.6| ) we have 

Ml = m]j + An, (B.7) 

where we have defined (we drop the subscript R from IIr) 

An = n(p2 = M^)-n(p2 = o). (b.s) 

Note that m^j defined in (|B.6| ) is renormalization scheme dependent, as V^ff is, while 
Mfj is not. In particular, in eq. (|B.7|) both m^ and AH depend on the renormalization 



scale /i in such a way that M'jj results to be scale independent (at least to 0{h)). In 
fact, from eq. (p.3|) and the definition 



IdlogZn 
^ 2 dlog/i ^ ' 

we easily obtain 

-^ = -27Mi (B.IO) 

which cancels the scale dependence of mjj (see eq.( ^7^) ). 

We now want to give the complete expression for the physical mass Mfj. The 
quantity All in the Landau gauge is given by the sum of the following terms 

An = AUft (top contribution) 
+ A]Iiv±iv^ + AUzozo 
+ An^^i^T + An^o^3 (gauge and Goldstone bosons contribution) 

+ An^i^T + An^3^3 

-|- AUhh (pure scalar bosons contribution). (B.ll) 
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In eq. (|B.11| ) we have taken into account only the contribution from the heaviest 
fermion, the top, and indicated by x^ a^id Xa the charged and the neutral Goldstone 
bosons, respectively. 

The complete expression for the different contributions to All calculated in the MS 
scheme is, for Mu < 2My/: 

i) Top contribution: 



An,i 



87r2 



-2Mt 



'Mf 
Ml, 



2 



H 



\ Ml 
log^ 

/i2 



'mT 

Ml, 



[B.12) 



ii) Gauge bosons and Goldstone bosons contribution: 



An^yiiyT + AUzoz" + AIIvi/i^T + An^o^3 
g'^M^ ( ^(1) ^ ^(2) , Ml (3) 






^TT" 



w 
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(B.13) 



where all the masses in the above expression have to be understood as the physical 
ones. The 11^ {i = 1, ....,9) functions read 
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where Z(x) is the function 




2Atan-i(l/A), 
Alog[(l + A)/(l 



if a; > 1/4 
A)] , if a; < 1/4 



[B.15) 



The terms containing the factor log{l — Ml / Ml ^) develop an imaginary part in the 
region M^z < Mu < 2My/ corresponding to the unphysical decays H — > W^x^ ^ Z^X^- 
However, this imaginary part, along with the whole factor log(l — Ml/Ml^), cancels 
in (B.13). In fact, using (B.14), (B.13) can be written as: 
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There is also an explicit imaginary part in (B.16) giving rise to 



Im(Anvi/±H/T + An^o^o 
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1287r M^ 

which will also cancel, as we will see. 

ill) Pure scalar bosons contribution: the contribution to AH coming from the pure 
scalar sector deserves more attention and we want to discuss it in more details. 

It is well-known that in the Landau gauge the Goldstone bosons x's do have a field 



dependent mass rn^{(f)) 



-m 



' + X(j)/2 which vanishes at the minimum of the potential 



V^ff((/)). As a consequence, the running mass niH presents an infrared logarithmic 
divergence when Goldstone bosons are included in the effective potential V^g ((/)). On 
the other hand, the physical mass Mh must be finite and gauge independent, so the 
divergent contribution coming from the Goldstone bosons to m?^ must be cancelled by 
an equal (and opposite in sign) contribution of the same excitations to All. To see 
it explicitly, one can imagine the Goldstone bosons to have a fictitious mass m-^ and 
calculate their contribution Am'jj to the running square mass mfj. It is not difficult to 
see that this contribution from the effective potential is 
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where the scale /i(t*) is defined in the text. 
An from the pure scalar sector reads 



On the other hand, the contribution to 
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where the last term comes from the Feynman diagrams involving Goldstone bosons. 
The explicit imaginary part in (B.19) would correspond to the unphysical decays H —>■ 
X^X^^ X^X^ and cancels against eq. (B.17). Expanding now the function Z (m'^/Mlj 
around m^ = one can easily show that the logarithmic divergence in Ami disappears 
and the final result for the pure scalar bosons contribution to AMI is finite and given 

by 

3 9^M'„ 



AMI 
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TT 



V^ 



Ml 



(B.20) 



1287r2 

Finally it is worth making a couple of comments. First, we have included the pure 
scalar bosons sector in AH for the sake of completeness, but now one is no longer 
allowed to compare the physical mass given in eq. (|B.7D with the running mass since 
the latter (see eq. (p.l8| )) is not well defined when Goldstone bosons are taken into 
account. Secondly, there is no scale dependence in the last expression ( [B.20|) (the scale 
n{t*) is fixed (see the text)), in agreement with the fact that the A-dependence of the 
anomalous dimension 7 of the Higgs field arises at two-loop. 
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Figure Captions 

Fig. 1 Plot of (l)mmit)/C,(t) as a function of fi{t) for rrit = 160 GeV and supersymmet- 
ric parameters: Ms = 1 TeV, Xt = 0, cos^ 2/3 = 1. Tlie solid (dashed) curve 
corresponds to the two loop (one-loop) approximation. 

Fig. 2 Plot of fi{t*) as a function of rrit for cos^ 2/5=1 (solid line) and cos^ 2/3 = 
(dashed line) and values of supersymmetric parameters: a) Ms = 1 TeV, Xt = 0; 
b) Ms = 10 TeV, Xt = 0; c) Ms = 1 TeV, X^ = QMl; and d) Ms = 10 TeV, 
X^ = 6M|. 

Fig. 3 Plot of mH,der(t)/mH(t) as a function of fi(t) for nit = 160 GeV and the super- 
symmetric parameters: Ms = 1 TeV, Xt = 0, cos^ 2/3 = 1. 

Fig. 4 The thin lines correspond to ttlh as a function of rrit in the two-loop (solid) 
and one-loop (dashed) approximation for supersymmetric parameters: Ms = 1 
TeV, Xt = 0, cos^ 2/? = 1. The thick solid line is the plot of Mh, (4.3), as a 
function of Mt, (4.1), in the two-loop approximation and for the same values of 
the supersymmetric parameters. 

Fig. 5 Plot of the physical, M^ (solid curve), and running, rriH (dashed curve) Higgs 
mass as a function of the scale /i(t) for supersymmetric parameters Ms = 1 TeV, 
Xt=0, cos^ 2/3 = 1. 

Fig. 6 Plot of Mh as a function of Mt for cos^ 2/3 = 1 (sohd line) and cos^ 2/3 = 
(dashed line) and values of supersymmetric parameters: a) Ms = 1 TeV, Xt = 0; 
b) Ms = 10 TeV, Xt = 0; c) M5 = 1 TeV, X^ = 6M|; and d) Ms = 10 TeV, 
X? = 6M|. 

Fig. 7 Plot of Mh as a function of tan/3 for Mt = 170 GeV and Ms = 1 TeV. The 
solid (dashed) curve corresponds to the Xf = 6M| {Xt = 0) case. 

Fig. 8 Plot of Mh as a function of /i(t*) for Mt = 170 GeV, Ms = I TeV, Xt = 0. 
Solid (dashed) lines are the two-loop (one-loop) approximation, thick (thin) lines 
correspond to cos^ 2/3 = 1 (cos^ 2/3 = 0). 

Fig. 9 Plot of Mh as a function of Mt for Ms = 1 TeV and tan/3 determined from 
the fixed point condition ( [S.l] ). The solid line corresponds to maximal mixing, 
Xf = 6M| and the dashed line to Xt = 0. 
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